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1.  Introduction 
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The  design  of  maintenance  policies  for  maintainable  and  repairable 
systems  makes  use  of  information  about  probability  distributions  of  com- 
ponent lifetimes  and  downtimes.  Since  such  designs  are  facilitated  if  these 
distributions  have  simple  analytic  forms,  downtime  distributions  have  been 
frequently  modeled  as  lognormal,  Weibull  or  Erlarg  in  form^  These  dis- 
tributions all  are  skewed  and  correspond  to  non-negative  random  variables 
like  downtimes.  Here  we  consider  one  more  family  of  distributions  which 
has  some  physical  motivation. 

Since  a downtime  interval  is  often  the  sum  of  subsidiary  intervals 

(for  failure  isolation,  component  removal,  repair,  reassembly,  alignment, 

etc.  ) it  seems  reasonable  to  think  of  the  downtime  x as  a sum  of  subsid- 

n 

iary  time  intervals 


x 

n 


r. 

1 


(1.  1) 


The  subscribt  on  x reminds  us  of  the  number  of  summands.  Several  dis- 
n 

tributions  are  possible  for  the  individial  r^,  but  here  we  consider  exponential 
distributions  which  are  the  simplest  and  which  are  widely  used  to  represent 
random  times  between  events. 

It  is  well  known  that  the  down-time  x will  have  an  Erlarg  distri- 

n 

bution  if  the  r^  are  independent  exponential  variables  with  identical  mean 
values.  Muth  [ 2 ] has  considered  the  approximation  of  Weibull  and  log- 
normal distributions  by  x^  in  which  the  r^  are  independent  exponential 
variables  but  with  possibly  diffe  rent  mean  values.  Here  we  further  gener- 
alize to  allow  de pendence  among  the  r^ a reasonable  situation  if  the 

vairables  represent  related  steps  in  a sequence  of  downtime  operations. 


1 


Section  2 discusses  several  forms  of  multivariate  distributions 
which  have  exponential  mariginal  distributions.  The  one  based  on  the  sums 
of  sequences  of  normal  variables  is  selected,  because  of  its  analytical 
simplicity,  for  use  in  the  Section  3 calculations  of  downtime  distributions. 

In  addition  to  giving  examples  of  these  cor  related-sums  distributions.  We 
show  that  introduction  of  dependence  among  r^,  (with  possibly  unequal 
means)  does  not  broaden  the  class  of  xn  distributions  over  that  which  re- 
sults from  independant  r^.  That  is,  the  sum  of  n dependent  exponential  vari- 
ables has  a distribution  identical  to  that  of  the  sum  of  n other  independent 
exponential  variables. 

Concluding  comments  refer  to  approximation  of  lognormal  variables 
by  sums  of  exponential  ones,  along  with  other  possible  extensions. 

2.  Multivariate  Exponential  Distributions 

It  is  well  known  that  a multivariate  distribution  is  not  uniquely 
specified  by  its  marginal  distribution^^.  For  example,  while  the  bivariate 
normal  distribution  has  nice  analytical  properties,  it  is  not  the  only  one 
with  normal  marginals.  This  multiplicity  will  be  demonstrated  for  the  case 
of  exponential  marginals  by  considering  a few  possible  bivariate  densities. 

Gumbel  [ 3]  considered  several  bivariate  exponential  densities. 

The  first  F.(rlf  r.,)  is  based  on  the  following  general  formula  for  com- 
bining marginal  distributions: 

F (x,  y)  = Fx(x)  Fy(y)  | i + a [ 1_FX(X)  1 l 1-Fy(y)  1 j (2.  1) 


When  applied  to  exponential  variables  with  unit  mean  values  this  produces 


**+1  ML 


' ■MfiyiWii 
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the  distribution  function 


"rl  'r2  ‘rl"r2 

Fjtrj,  r2)  = ( 1-e  *)  (1-e  C)  (1  + a e C) 


(2.2) 


rl’  r2  - ° 


la  I < 1 


and  the  density  function 


-ri-r2  r -ri  -r2 

f1(r1,r2)  = e '[l  + a(2e  -1)  (2  e J-l) 


(2.3) 


In  this  model,  a - 0 corresponds  to  independence  of  r j and  r2,  and  it  can 

be  shown  that  the  correlation  coefficient  p is 

rl’r2 


rr  r2 


= a/4 


(2.4) 


with  its  magnitude  limited  to  be  less  than  1 / 4. 

Another  model  of  Gumbel's  is  defined  by  the  distribution  function 


. , -rl  "2.  -rl-r2-0rlr2 

F2(r  j,  r2)  = 1-e  -e  + e 


(2.5) 


rl  , r2  > 0 


0 < 0 < 1 


and  the  corresponding  density 


ri  r2  e ri  rz[(i  + e r^d+e  r2)  -e]  (2.< 


f2(rl'r2,=  e 


nun  *<  ;* 
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Here  0 = 0 corresponds  to  independence,  and 


p Q- 1 + l/0  „ . 1 , , 

r j , r ^ = - 0 e ' Ei  (0  ) - 1 

( Ei  is  Jn  e* jaenerttia  I 

Marshall  and  Olkin  [ 4 ] introduced  the  bivariate  exponential  de- 
fined by 


F (rl’  r2)=P  fRl>rl'  R2  >r2] 


"X  lr  1 " X 2 r2  " X 12  max^r  r2^ 


(2.7) 


F3(r1.r2)=  1 - e'|Xl  + ^2)rl  -r2<X2  + ^12> 


“rl  X 1 " r2  K 2 " K 12  max(r  1'  r2) 

+ e 


(2.8) 


Here  the  correlation  coefficient  is 


Prl->r2  = X12/(X1  +X2+X12) 


(2.9) 


which  ranges  between  zero  and  one.  This  distribution  can  be  thought  of 
as  a result  of  fatal  shocks  occuring  from  three  independent  Poisson  sources 
with  rates  X.  j , X.  ^ and  X 12*  Component  1 with  lifetime  rj  is  killed  by 
events  of  the  first  or  third  variety,  and  r^  is  determined  by  events  of  the 
second  or  third  type.  Barlow  and  Prochan  [ 5 ] point  out  that  (2.8)  is  the 
unique  exponential  bivariate  with  the  zero  - memory  property:  The  joint 
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survival  probability  of  a pair  of  components  each  of  age  t is  the  same  for 
all  t (e.g.  the  same  as  if  both  were  new).  This  zero  memory  property  is 
quite  desirable  when  modeling  joint  lifetimes  of  components  in  a system. 

In  the  present  context  of  r^  representing  durations  of  a sequence  of  related 
events,  this  special  property  seems  inessential. 

Kibble  [6  ] considered  a bivariate  exponential  density  of  the  form 


F4(rl’r2;P)=  —4— 27 
4a  ( 1 - p ) 


'/.-Z 1 n — ' 

^Vp2)  Ie  JLLl/2  (2.1( 

L,J  (1-p  )J 


in  which  I0  is  a modified  Bessel  function  . 

We  will  derive  this  in  an  alternate,  more  natural  way  for  use  be- 
low. The  density  in  (2.  10)  applies  when  we  view  r j and  r ^ as  being  gene- 
rated from  correlated  normal  variables.  It  is  well  known  that  if  Wj  and 
Zj  are  independent,  zero  mean,  equal  variance  (cr ^ ) normal  variables  then 
r 1 defined  as 


2 x 2 
rl"wl  +Z1 


(2.  11) 


has  an  exponential  distribution  with  mean 


E (r.)  = 2 a‘ 


(2.  12) 


Now  if  (wj,  w^)  and  (Zj,  z^)  are  two  independent  pairs  of  normal  variables. 


but  with 


cov  (w.,  w.)  * cov  (z.,  z.) 
1 J x J 


i ■ 1,2 


j * 1.  2 


(2.  13) 
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then  Tj  and  r 2 defined  by 


2 . 2 . , , 
r.  = w.  + z.  1=1,2 

ill  * 


(2.  14) 


will  be  dependent  exponential  variables. 

Equation  (2,  10)  can  be  derived  from  this  reasoning  when  all  four 

2 

variables  ( Wj,  w^,  z^,  z^)  have  equal  variances  a . In  that  case 


_ 1 [ Wi2-2PW1W2+W22  + Zl2-2PZ1Z2+Z22 

f(wl’  W2»  Zl'  Z2  ” .2  2.  4 exp  " 2.  2 

(2ir ) (1  - p ) cr  L 2(l-p  ) <r 


(2.  15) 


pW2,  Zj,z2)  dw  j , dw2»dz  y dz. 


2 2 . 

1W1  + Z1  < rl 


2 2 

w2  + z2  < rz 


(2.  16) 


Introduction  of  polar  coordinates  in  Wj,Zj  and  w2#  z2  planes 


j = Yj  cos  w2  = Y2  cos  ®2 


;l  s Yj  sin  9j  z2=v2sin02 


Ur 
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reduce  the  F integral  to  the  form 


i ^*2  2ir  2ir 

F4(ri»  r2)  = ' 2 2 4“  Iff  I 

4ir  ( 1 - p j o-  4 Yl=0  y2=0  V0  e2=l 


-[\i  + Y22  - 2p  Yj  Y2  cos  (6^62)  / 2 <r2  ( l - p 2 ) J 


(2.  17) 


* Yi  y2  dYl  dy2  de}  de2 


Substitution  of  ^ = ” ®2  f°r  the  0j  integration  produces  a periodic  in- 

tegrand, independent  of  0£.  Thus  the  6 - integrals  become 


_ 2tt  2tt 

I J = 2lr  / 


2 2 

2tt  P Y1  Y2  cos  V ! * (‘"P  ) 


V°  V0 


2 |p|  Y1  Y2 
4ir  Irt  2,,  2. 

0 0-  (l-p  ) 


(2. 18) 


Finally  substitution  of  (2,  18)  into  (2.  17)  and  differentiations  with 

respect  to  r and  r..  produces  (2.  10).  It  turns  out  that  p is  just  the 

2 ri’r2 
square  p of  the  correlation  of  the  underlying  normal  variables,  with  a 

range  between  zero  and  one. 

The  Bessel  function  form  for  this  bivariate  exponential  based  on 
normal  variates  may  not  appear  to  be  very  felicitous.  However,  this  kind 
of  distribution  will  be  convenient  when  we  concentrate  in  the  next  section 
on  the  sum  of  dependent  exponential  variables,  as  in  (1.  1). 

For  the  n- dimensional  version  of  this  class  of  distributions,  we 
consider  zero  mean,  normal  n-vectors  w and  z each  with  the  same  covari~ 


FI 


ance  matrix 


El  w w1]  = E[  z zq  = r 


(2.19) 


In  this  way,  for  each  i w^  and  z^  will  have  the  same  variances  so  the  sum 
of  their  squares  will  be  an  exponential  random  variable  r^..  We  do  allow 
r.  and  r^  to  have  unequal  means,  contrary  to  the  special  case  in  (2.  10). 

3.  Examples  and  Properties  of  Proposed  Family  of  Distributions 

The  previous  sections  introduced  the  multivariate  exponential  dis- 
tribution based  on  a pair  of  multivariate  normal  distrubtions.  General  ex- 
pressions for  this  distribution  were  not  presented,  but  we  are  only  interested 
in  the  distribution  of  the  sum  of  the  dependent  exponential  variables. 

Some  general  ideas  about  such  distributions  are  revealed  in  the 
following  simple  examples.  We  consider 


x2  = ri  + rZ 


(3.  1) 


the  sum  of  two  equal  mean  exponentials  which  have  the  joint  density  in  (2.  10). 
As  we  will  see  in  a more  general  context  below,  the  characteristic  function 
of  X£  is 

*,<■>«  El  e"8*2] 


4(l  - p 2)  s2  + 4 s + 1 


(3.2) 


when  E[  r.  ] = E [ ] = 2 a * 2.  Table  I shows  several  examples  of 

f as  p varies  between  0 and  1.  (Recall  the  p is  the  correlation  coefficent 
x2 

c 2 
of  the  underlying  normal  variables  and  p is  the  correlation  coefficient  of 


.wfJHBgai  .# 


r j and  r^).  In  each  case  £[  ] * 4,  and  the  peak  of  the  density  moves  left 

toward  the  limiting  (p  * 1)  exponential  case  as  p increases. 

Table  I reveals  another  interesting  fact.  This  is  that  for  each 

value  of  p relating  r ^ and  r^  we  can  always  find  two  inde pendent  exponential 

variables  r^  and  r with  possibly  different  mean  values,  such  that  x^  = 

r'j  + r'^  has  the  same  distribution  as  x2  = rj  + r2* 

Another  interesting  comparison  is  that  of  x^,  . . . based  on 

addition  of  more  and  more  dependent  exponential  random  variables.  The 

natural  extension  of  the  bivariate  example  in  (3.  2)  is  to  consider  all  under- 

2 

lying  normal  variables  to  have  variances  of  c r , and  with  the  corresponding 
w and_z  vectors  each  having  the  covariance  matrix 


(3.  3) 


This  implies  that  the  r^  have  the  correlation  coefficient 


P * P 

r r.  r.  r 

i J 


2 |i  - j f 


(3.4) 


if 


>V-- 


Ik-. 


If  we  further  consider  mormalized  variables 


x = 
n 


1 


2n  <r 


r. 

1 


(3.5) 


then 


E[  x ] = 1 
L n J 


(3.6) 


Calculation  of  other  moments  for  these  x is  straightforward, 

n 

though  tedious,  from  the  moment  generating  properties  of  the  underlying 

normal  distributions.  Moment  expressions  for  the  special  example  using 

(3.  3),  given  in  the  appendix,  show  the  following  limitations  on  the  possible 

2 

values  of  variance  <r  , Skewness  t,  , and  Kurtosis  v • 

n n 11 


_1_  _<  or  < 1 

n 


2 a < f,  < 2 
n — n 


(3.7) 


6 <r  < v 
n — n 


< 6 


Figure  1 shows  the  effects  on  the  density  function  shape  of  increas- 
ing the  number  of  equal  mean  exponential  variable  which  are  summed. 

Those  curves  are  for  x ^ x ^ each  with  unit  mean  and  variance  of 
cr^2  * 0-2  9 - 0.  65  364583  (corresponding  respectively  to  p = 0.55433895  and 
p = 3/4).  The  equations  for  those  densities  were  found  to  be 


f (t>=  1.  8039505  | e'1*  2867206t  „e  - 4.  4877! 56t> 
x2 


f (t)  = 0.  99l59025e  -l6,  3l0l64t  + 

x3 


(3.8) 


• v.-.-  ■*. 
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+ 1.6750764 e 


■ 1 . 26 1 2644t 


- 2. 666667  e 


■ 6.  8571 428t 


The  similarity  of  the  two  curves  in  the  figure  attests  to  the  need  for  great 
precision  in  expressing  the  densities  and  in  computing  typical  points. 

Figure  1 also  shows  the  lognormal  density  with  the  same  mean 
and  variance  as  5.^  and  x 

The  relation  between  families  of  distributions  which  was  suggested 

by  the  results  in  Table  I will  now  be  generalized  in  the  form  of  a Theorem. 

We  define  r , r.,,  . . . r as  sums 
12  n 


2 . 2 
r.  = w . + z . 

ill 


(3.9) 


in  terms  of  two  independent  normal  n-vectors  w and  z , each  with  zero  mean 
values  and  covariance  matrix  T.  As  shown  above,  the  r^  are  exponentially 
distributed  with  mean 


E[  rj  ] = 2 Yii 


(3.10) 


and  correlation  coefficients 


'r,r,  * (V‘ 


(3.  11) 


..  u * 
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= (^)  'n  / exp  [-  j L'  R_1  l]  d L (3.i7> 

where  we  have  defined  the  matrix 

R"1  = I + 2 s M'  M (3.  18) 

The  integral  in  (3.  17)  is  of  the  form  of  a normal  density  integrated  over  all 
values,  except,  for  a scale  factor. 


I 

s 


(3.19) 


and  the  desired  characteristic  function  for  the  sum  of  exponentical  variables 
is 


(3.20) 


(several  examples  are  found  in  the  appendix). 

The  roots  of  q(s)  are  negative  reciprocals  of  the  eignvalues  of  the 
real  symmetric  matrix  2 M1  M,  so  they  are  real  numbers.  If  the  r^  are 
independent  then  M is  diagional 


and  so  is 


Mj  = diag  ( ccj , • • • orn) 


(3.21) 
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2 

qj(s)  = it  (1  + 2 a.  s)  (3. 23) 

i=  1 

The  proof  is  completed  by  comparing  the  polynomials  q(s)  in  (3.  20), 

for  a general  M,  to  the  qj(s)  in  (3.23)  for  independent  r^.  The  mean  values 
2 

of  the  latter  (2a.)  can  be  chosen  to  make  the  roots  of  qj(s)  match  any  n 
( necessarily  real)  roots  of  q(s).  The  resulting  polynormals  and  character- 
istic functions  will  be  identical  because  the  ®(0)  = 1 property  removes  any 
scale  factor  ambiguity.  This  completes  the  proof. 

In  general,  the  n-independent  exponential  variables  whose  sum  is 
indistinguishable  from  the  sum  of  n-  correlated  ones  will  have  different 
mean  values  from  those  of  the  correlated  variables. 

Once  this  structure  has  been  established  for  sums  of  correlated 
exponential  variables,  previous  results  for  sums  of  independent  variables 
(e.g.  those  of  Muth  ) are  directly  applicable.  However,  it  appears  that 
when  equal  mean  variables  have  their  mean  and  correlation  p in  (3.  3)  ad- 
justed so  xr  matches  the  lognormal  mean  and  variance,  then  f tends 
to  lognormal  for  large  n.  This  is  in  contrast  to  Muth's  summing  of  inde- 
pendent variables  in  which  he  had  to  search  for  the  proper  individual  mean 
values  which  made  the  sum  approximately  lognormal. 

Other  interesting  properties  follow  from  knowledge  of  the  general 
form  of  the  solution.  For  example,  if  all  roots  of  (11)  are  distinct  then 

n - X.  .t 

fx(t)=  S a.e  1 (3.24) 

n i=  1 

(The  roots  must  be  negative  since  f(t)  must  have  finite  area).  One  condition 
on  the  a.  is  that 


I 


<tt 

f (t)  = 

X A 

n 


(3.2'.; 


* 


application  of  the  initial  value  Theorem  of  Laplace  transforms  to  (3.20) 
shows  further  that 


f 00  (0)  = 0 k * 0,  1,  ..  (n-2)  (3.26) 

n 

n 

= S (-X.)*  a. 

1 1 1 

Equations  (3.25)  and  (3.26)  provide  n linear  equations  which  can  be  solved 

for  the  n coefficients  a^  in  (3.  24),  if  the  roots  X^  are  known. 

4.  Comparisons  with  Lognormal  Distribution 

A lognormal  random  variables  x is  defined  by  saying  that  its  natural 

2 

logarithm  has  a normal  distribution  with  mean  r/  and  variance  <r  . Thus, 
the  density  of  x is 

-1  - (In  x - n)2/  2 o-2 

f (x)  = (xa  V 2-rr  ) e ; x _>  0 (4.1) 

and  its  mean  and  variance  can  be  expressed  as 

Trf  i p + cr2/  2 

E[xJ  = e ' 

Var  [x]=e2  n + ^ [ e ^ -\  ] 

- E2  l x]  [ e * - l] 

Figure  1 shows  a typical  lognormal  density,  which  is  necessarily  defined 
for  positive  x and  is  unimodal. 

Lognormal  models  have  been  used  frequently  to  characterize  down- 
time data,  probably  because  of  the  ease  in  doing  so  with  the  aid  of  normal 
probability  paper.  Many  other  simple  unimodal  distributions  can  often  be 
used  with  equally  good  effect  with  a finite  amount  of  data. 


(4.2) 

(4.3) 
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One  special  property  of  the  lognormal  is  its  behavior  for  large 
arguments.  One  way  to  describe  this  behavior  is  in  terms  of  the  conditional 
mean  exceedance 

00  _ 

CME  = E[  t - t | t>  t ] s J F (t)  dt  / F(t  ) (4.  4) 

T 

If  t were  a lifetime  random  variable,  then  CMF  would  be  the  expected  future 
life  as  seen  at  time  t.  CME^  is  a constant  for  an  exponential  random  vari- 
able, but  is  an  increasing  function  of  t,  when  t is  large,  for  a lognormal 
random  variable.  Thus,  the  latter  is  said  to  have  a heavy  tail. 

It  is  clear  that  the  tail  of  an  x as  defined  in  the  previous  sections 

n 

will  not  be  heavy,  since  each  density  there  has  the  form  of  a sum  of  decay- 
ing exponentials.  The  tail  of  such  a density  will  be  dominated  by  the  slow- 
est of  these  exponentials,  so  it  will  be  neither  heavy  nor  light. 

6 . Conclusions 

Several  multivariate  exponential  distributions  were  examined  for 
use  in  modeling  constituent  time  intervals  whose  sum  represents  a system's 
downtime.  The  distribution  based  on  sums  of  squares  of  pairs  of  normal 
variables  was  chosen  for  its  analytic  simplicity. 

We  showed  that  variation  in  the  correlation  and  means  of  such 
variables  does  not  produce  any  distributions  for  the  sum  which  could  not 
have  been  realized  by  a sum  of  inde pendent  exponential  variables.  Such 
sums  can  be  found  to  well  approximate  the  frequently  used  lognormal  down- 
time distribution,  except  for  the  tail  behavior.  It  is  not  clear  that  the  heavy 
tailed  behavior  of  the  lognormal  is  truly  more  appropriate  for  down  time 
modeling,  A lognormal  model  could  be  justified  for  a variable  resulting  from 
a produc  t of  independent  variables  (by  applying  the  Central  Limit  Theorem 


to  its  logarithm),  but  there  seems  to  be  no  physical  justification  for  such 


a downtime  model* 


The  multivariate  exponential  model  studied  here  may  also  be  useful 
in  generalizing  maintenance  and  replacement  strategy  modeling  by  introduc- 
ing dependent  intervals  into  problems  where  independent  variables  have 
been  used  previously.  It  is  interesting,  in  this  context,  to  note  that  when 
a development  paralleling  that  of  (2.  10)  is  applied  to  sets  of  three  correlated 
normal  variables  with 


E [ w.  w.  ] = O'.  <r.  p ^ “ j I 
1 i J J i J 

then  the  trivairate  exponential  density  becomes 


/ 2 

1 2 

Vrir2p 

T 

Jr2r 3 P 

2. 

I0 

/.  2. 

Q» 

1 

'N) 

b 

L 
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APPENDIX 


Moments  for  Sums  of  Dependent  Exponential  Random  Variables: 
Case  of  Equal  Means  and  Exponential  Correlation  Coefficient 


This  appendix  gives  further  details  for  the  distributions  introduced 
in  Section  3 and  defined  by  (3.  3)  with  normalization  (3.5).  We  consider 
cases  for  n x 2,  3,  and  4 using  (3.  20)  get  the  characteristic  functions 


* (s)  * [ 0.  25(1- p2)  s2  + s + l]  " 1 


*3(s)  *[(8-l6p2  + 8p4)(s/6)3+(l2-8p2-  4 p 4)  (s/£)2  + s + 1 ] 


S (s)  = [(l6--48p2  + 48p  4 - l6p6)  (s/8)4  + (32  - 48p2  + l6p6)  (s/8)' 


+ (24-  12p  2 -8p4-4p6)(s/8)2+s  + l]  _I 


Moments  of  interest  are  defined  as  follows: 


m (n)  = E f x K I 
K 1 n 1 


mean  = u = m j (n) 

Variance  = <r  2=m,(n)  - m2  (n) 
n c.  l 


Shewness  = L » E 
n 


Kurtosis  = v = E 
n 


Xn  " mi  (n)  H 
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- 3 3 

= it  [ m^(n)  - 3 m^n)  m^(n)  + 2 m ^ (n)  ] 


<t~^  l rr^(a)  - 4 rrij  (n)  m^(n)  + 6 m^njm^n) 


- 3 m j (n)  ] - 3 


Formulas  for  these  moments  as  functions  of  p are  given  below: 


n = 2 


= 1 


2 1,1 
* * 2 + 2 P 


£ = 


,1  4.  1 n2l3/Z 

(2  + 2 P ’ 


(0.5  + 1.5  p6) 


(0.5  + 0.5  p ) 


2 ^ (1*54  6p  2 + 1. 5 p 4 ) 


n=  3 


= l 


2i  4 2 2 4 

<r  = 3-  + ? P + 9 P 


; = 


1,4  22  4.  3/2 

(j  + 9 P + o P ) 


(0.  222222  + 0.888888  p2  + 0.888888  p4) 


v = (0.55555  + 2.  66  565  p2  + 3.  501673  p 4 + 1.  185 1 85  p 


(<r  ) 


+ 0.296296  p8)  - 3 


0.  25  + 0.  375  p + 0.  25  p + 0.  1 25  p 


VsTT*0-  125  + 0 


~TT  (°*  28125  + 
[<r  ) 


0. 7529297 


5625  p2  + 0.  75  p4  + 0.  375  p6) 


1.40625  p2  + 2.8125  p4  - 0.  28125  p6  + 
p8  + 0.  375  p 10  + 0.  09375  p12)  - 3 


— - . .. 


Table  I : Densities  for  Sums  of  Exponential  Variables 

(a)  Equal  mean  and  Correlated 


(bj  Independent  with  Unequal  means 


